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We analyze the crosstalk error mechanism in measurement of two capacitively coupled supercon- 
ducting flux-biased phase qubits. The damped oscillations of the superconducting phase after the 
measurement of the first qubit may significantly excite the second qubit, leading to its measurement 
error. The first qubit, which is highly excited after the measurement, is described classically. The 
second qubit is treated both classically and quantum-mechanically. The results of the analysis are 
used to find the upper limit for the coupling capacitance (thus limiting the frequency of two-qubit 
operations) for a given tolerable value of the measurement error probability. 
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I. INTRODUCTION 

Superconducting Josephson-j unction circuits, includ- 
ing phased flux)S and charge'^ qubits, have attracted a 
significant interest as promising devices for quantum in- 
formation processingi^ In this paper we consider flux- 
biased phase qubitSfS*Sil which have been introduced rel- 
atively recently and have a clear advantage over the 
current-biased phase qubits. While the schematic of a 
flux-biased phase qubit (Fig. ^ may be very similar to a 
flux qubit (in the simplest case, a superconducting loop 
interrupted by one Josephson junction), an important 
difference is that in the phase qubit the logic states are 
represented by two lowest levels in one well of the cor- 
responding potential profile, while for the flux qubit the 
levels in two neighboring wells are used. An imaginary- 
swap quantum gate, which together with single-qubit ro- 
tations forms a universal set of quantum gates^^ has been 
realized with flux-biased phase qubits in Ref . '!J. 

Simultaneous measurement of all the qubits is an im- 
portant stage in quantum information processing, rep- 
resenting one of numerous challenges on the way to a 
scalable quantum computer. Recently a scheme allowing 
for fast simultaneous measurement of two coupled flux- 
biased phase qubits has been implemented.^ According to 
this scheme, a measurement is performed by lowering the 
barrier between the shallow ("left") potential well used 
for qubit states and a much deeper "right" well (Fig. 
P), so that a qubit in the upper state |1) switches by 
tunneling to the right-hand well with probability close 
to one, whereas a qubit in the lower state |0) remains 
intact. Since the qubit coupling is fixed in the sim- 
ple design realized experimentally so far, the measure- 
ment is complicated by a crosstalkiS The mechanism of 
the crosstalk is the following. Suppose that one of the 
qubits is measured in state |1), which means tunneling 
to the right well. The tunneling will then be followed by 
the dissipative evolution (damped oscillations), eventu- 
ally reaching the ground state of the right well. These 
oscillations will obviously perturb the state of the sec- 
ond, capacitively coupled qubit, especially because in the 




FIG. 1: The circuit schematic of a flux-biased phase qubit 
and the corresponding potential profile (as a function of the 
phase difference 5 across the Josephson junction). During 
the measurement the state |1) escapes from the "left" well 
through the barrier, which is followed by oscillations in the 
"right" well. This dissipative evolution leads to the two-qubit 
crosstalk. 



experiment the oscillation frequencies in both wells are 
somewhat close to each other, and both qubits are practi- 
cally identical. Therefore, if the measurement of the sec- 
ond qubit is made after the dissipative evolution of the 
first qubit, the measurement result is likely to be wrong 
(it will most often be |1) because of the second qubit 
excitation^). To avoid this problem, the measurement 
of both qubits should be done almost simultaneously^^ 
within the timescale shorter than development of the 
crosstalk mechanism. However, this is not a complete 
solution of the problem because the excitation of the sec- 
ond qubit due to crosstalk may be sufficient to switch the 
qubit from the left well to the right well, even if this ex- 
citation happens a little after the second qubit measure- 
ment, that will also lead to the measurement error. We 
emphasize that this crosstalk mechanism is due to a fixed 
coupling between the qubits. It seems possible to realize 
schemes with adjustable coupling in future, which will 
eliminate much of the measurement crosstalk discussed 
here. Still, it is of interest to analyze the fixed-coupling 
case, since that is the simplest scheme and the only one 
realized experimentally so far. 

The present paper is devoted to the theoretical study 
of the measurement crosstalk between two capacitively 
coupled flux-biased phase qubits. In Sec. ^ we study 
some basic properties of phase qubits. In particular, the 
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simplified cubic potential is discussed and the Hamilto- 
nian for two capacitively coupled qubits is derived. The 
measurement crosstalk is studied in Sees. IIIIIVI For 
definiteness, we assume that the first qubit is switched 
(state |1)), whereas the second qubit is initially in the 
state |0). The dynamics of the first qubit, which after the 
switching performs damped oscillations in the deep well, 
is analyzed classically in Sec. lIIII Such an approximation 
drastically simplifies the problem and, at the same time, 
is quite accurate, since for the experimental parameters^ 
used here the first qubit is highly excited after measure- 
ment, with typical quantum number over 10^. The sec- 
ond qubit in this paper is treated both classically and 
quantum-mechanically. The classical treatment ("Sec. lIVII 
involves two approaches: the harmonic-oscillator modelf^ 
which allows for an analytical treatment, and a numeri- 
cal solution for the exact Hamiltonian. In the quantum 
approach (Scc.0), an efficient method of numerical solu- 
tion of the Schrodinger equation, based on using a sub- 
set of the eigenstates of the unperturbed Hamiltonian, 
is developed. The results of the quantum approach are 
somewhat similar to those of the classical treatment. Our 
quantum treatment, in contrast to the classical one, does 
not allow for dissipation in the second qubit, but an in- 
sight into possible effects of dissipation can be obtained 
by the comparison of the two approaches. The conclu- 
sions following from the present work are summarized in 
Sec. |VD 



II. FLUX-BIASED PHASE QUBITS 

Before the discussion of the basic equations for coupled 
qubits in Sec. Ill CI we review the properties of one qubit 
in Sees. ITTH and ITlBl 



A. Qubit potential 

A flux-biased phase qubit schematic^ coincides with 
that of the basic rf SQUII^ (Fig.[TJ. Neglecting dissi- 
pation, it can be describedi^ as a fictitious mechanical 
system with the Hamiltonian 



H 



2m 



U{S), 



(1) 



where i5 is the Josephson-junction phase difference, p — 
mS is the corresponding momentum, m = (<l>o/27r)^C 
is the effective mass determined by the capacitance C, 
$0 = h/{2e) is the flux quantum, e is the electron charge, 
and U{S) is the potential energy (shown schematically in 
Fig.UJ) 



UiS) = Ej 



2\ 



COSiJ 



(2) 



Here Ej = $o^o/27r is the Josephson energy, A = 
27r/oL/$o is the dimensionless inductance, 4> = 27r<f>/(f>o 




FIG. 2: The graphical solution of Eq. @; see text. 



is the dimensionless external magnetic flux, /q is the crit- 
ical current, and L is the inductance. 

In this and the next subsections we review the basic 
propertiesiS (see also Ref. 0) of the potential energy . 
Since the potential is invariant with respect to simulta- 
neous change of 6 and </> by 27rn where n is integer, we 
limit ourselves by the range < ^ < 27r for the external 
flux. Notice that the potential is symmetric for cf) = rtTr 
($ = n$o/2). 

The maxima and minima of the potential ^ satisfy 
the equation 



(J-0)/A 



(3) 



which can be solved graphically. Figure[21shows the r.h.s. 
of Eq. (PJ (curve 1) and the l.h.s. (straight lines 2-5) for 
several values of the parameters. This equation can have 
only one solution when the straight-line slope is greater 
than 1 (as for line 2 in Fig. O ; therefore in the case A < 1 
the potential (0) has only one well for any external flux 

For A > 1 the potential may have more than one well. 
There will be at most two wells if the slope of the l.h.s. 
of Eq. lO is greater than the slope 1/Ai of the line 3 in 
Fig. 121 which is tangent to — sin (5 at two points ±5' (the 
line passes through the origin because of the symmetry). 
This condition yields the equation tan 5' — 5' with the 
least positive root S[ = 4.493, which corresponds to Ai = 
-l/cos(5i = 4.603. Thus, for A in the intervaliS 



1< A < 4.603 



(4) 



the potential has one or two wells, depending on cj). In 
particular, for the experimental parameters of Ref. 9 used 
in this paper [see Eq. (|12|l below] one obtains A = 3.72 
(as for lines 4 and 5 in Fig. |2Jl, which satisfies condition 
Similarly, one can show that the potential will have 
n or n -|- 1 wells (depending on (jj) if 



A„_i < A < A„ {n> 1), 



(5) 



where Aq = 1, A„ = 1/| cos(5^|, and 6'^ is the nth (in the 
increasing order) positive root of the equation tan S' = S' . 
In particular, A2 = 7.790, A3 10.95, A4 = 14.10, and 
A„ {n+ l/2)7r for n > 1. 
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The condition for a two-well potential in the case Q 
can be found by considering the transition between the 
one- well and two- well cases, which is illustrated by lines 
4 and 5 in Fig. 13 These lines are tangent to — sin 6 at 
the points Sc (critical fluxes) which correspond to the 
inflection points of the potential |5J| (when a well disap- 
pears, the corresponding maximum and minimum of the 
potential merge, so that both the first and second deriva- 
tives are zero at this point). Solving the equation for the 
inflection points 



cos 6c = — 1 / A 



(6) 



(which does not depend on the external flux), we get two 
solutions in the interval (0, 27r): 



5c = tt/Z + arcsin(l/A), 5'^ = 2tt-5c 



(7) 



Inserting these results into Eq. lO, we finally obtain the 
condition for a two-well potential: 

0^<0<0c, (8) 

4>c = 7r/2 + - 1 + arcsin(l/A), 0c = 27r - 0c (9) 

(it is easy to show that ^J, < 0c for A > 1). 

For definiteness we consider the case tt < < 0c, 
in which the right well is deeper than the left well. The 
depth AUi of the left well (i.e., the difference between the 
potential maximum and minimum) can be characterized 
by the crude estimate of the number of discrete levels in 
the well 



N, 



AUi 

TlLOl 



(10) 



where is the "plasma" frequency (classical oscillation 
frequency near the well bottom) for the left well (simi- 
larly, u!r denotes the plasma frequency for the right well): 



^^Lr = ^JEj{1/\ + coaSi^r)/m; (11) 

here di^r corresponds to the well minimum and obeys Eq. 

Notice that Ni^^ is not necessarily integer and there 
is no simple relation between iV; and exact number of 
discrete levels in the well because of significant anhar- 
monicity of the potential. 

For numerical calculations presented in this paper we 
will use the following values of the parameters from the 
experiment of Ref. 

C = 700 fF, L = 0.72 nH, Iq = 1.7/^A. (12) 

Figure |31 shows the qubit potential U (6) for Ni — 10 
(corresponding to = 4.842), iV, = 5 (0 = 5.089), and 
Ni = 1.355 (0 = 5.308); the last value corresponds to 
the bias during the measurement pulse (see below) . The 
qubit levels |0) and |1) are, respectively, the ground and 
the first excited levels in the left well. 
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FIG. 3: The qubit potential U{S) [Eq. @] for Ni = 10 (dotted 
line), Ni — 5 (solid line) and A'^; — 1.355 (dashed line). 



B. Cubic potential 

When one of the wells is very shallow, it can be approx- 
imated by a cubic potential. Assuming < 0c — <C 0c 
(shallow left well) we can approximate cos 6 in the vicin- 
ity of Sc as 

cos S = cos Sc cos y — sin Sc sin y 

«-A-i(l-yV2)-yT^A^(y-j/V6), (13) 

where y = 5 — Sc- Then the potential (|2Jl can 
be approximated by the cubic polynomial Uc{S) = 
Vl — A^^ Ej{ey/2 — y^ /Q), where a constant is neglected 
and 



= 2(0c-0)//a^ 



1. 



(14) 



The minimum and maximum of this potential are at yi 
and j/max = \/e, respectively, i.e., at 



6l= 5c- \/e, (5max = Sc + 



(15) 



Shifting the axis as x = y + ^/e = S — Sc + \ft and again 
neglecting a constant, the potential can be rewritten as 



Uc{S) = yJ\-\-^Ej{^tx'l2 - xV6). (16) 
In this approximation^*^ the left well parameters are 



Ni = ^e'/\l-X-'y/*./^j. (17) 

The validity condition for the cubic-potential description 
()16|l is e <C 1, which is well satisfied for the left well with 
qubit parameters considered in the present paper. 



C. Two capacitively coupled qubits 

Let us consider two capacitively coupled flux-biased 
phase qubits^ (Fig.01). The current balance for this cir- 
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FIG. 4: The circuit schematic for two capacitively coupled 
flux-biased phase qubits. 



cult yields the equations 



C'Ri 



sin (5,- 



(18) 



where qubits are numbered by i, j — 1,2 (i ^ j); Cj' = 
Ci + Cx', the quantities 5i, /oi, Ci, and Li are respectively 
the Josephson-j unction phase difference, critical current, 
capacitance, and inductance for the ith qubit, Cx is the 
coupling capacitance, and (j)i — 27r$i/$o is the dimcn- 
sionless external magnetic flux. Dissipation in Josephson 
junctions is described here using the resistively shunted 
junction (RSJ) model^ by introducing resistances Ri 
and i?2 into the circuit (Fig. For each junction the 
dissipation can be characterized by the energy relaxation 
time 



Ti = i?iCi, Tl = i?2C'2- 



(19) 



In the absence of dissipation {Ri = oo) Eq. (|18|l can 
be written in the form of Lagrange's equationsi^ 



d^dL 



dL_ 



0, 



with the Lagrangian 

L = K-Ui{5,)~U2{S2). 
Here the ith qubit potential energy is [cf. Eq. ||2Jl] 

■(<5.-0,)' 



Ji 



2A,; 



COS Si 



(20) 



(21) 



(22) 



where Eji — $o^oi/27r, — 2nIoiLi/^0: and the kinetic 
energy K is 



K = 



m[Sl 



mxSiS2 



(23) 



where m- = (<i>o/27r)2C| and rrix = {^o/'^t^)^Cx are the 
normalized capacitances. Thus, the problem of two cou- 
pled qubits is equivalent to the motion of a fictitious par- 
ticle in the two-dimensional space {Si,S2)- From Eqs. 
(|21|l and (|23|l one can obtain the generalized momenta 
Pi = dL I dSi in the form 



Vi 



m'.Si 



rrixSj. 



(24) 



It is easy to see that pi is the total (node) charge on the 
ith junction and the adjacent coupling capacitor multi- 
plied by h/2ew^^ Dissipation can be accounted for by 
the addition^ of the friction force Fi into r.h.s. of Eq. 
(1201) ■ This yields Eq. ^ if 



F., = - 



$2 



Equations (|24|l can be inverted, yielding 



4 ^ Pt/m'l +pj/m'^, 



(25) 



(26) 



where 



= mi-\- (mj^ + iTT-x^)^^ and = nii + TO2 + 
mim2/mx- Inserting Eq. 12t)|) into (|23|1 . one obtains the 
Hamiltonian" H = K + Ui + U2 ro- the formic 



H = 



2 
Pi 

2m'/ 



2 

P2 ^PlPl^Ui{S^) + U2{S2). (27) 



Notice that the Hamiltonian H27|) can also be derived in 
a direct way (without Lagrangian language) using the 
facljiii that the node charge (muliplied by h/2e) pi is the 
conjugated variable to the phase Si and expressing the 
combined electrostatic energy of capacitors Ci, C2 and 
Cx in the form of three first terms of Eq. (|27|l . 

Instead of Lagrange's equations (j^Ol, one can use 
Hamilton's equationsii with respect to Si and pi, viz., 
Eq. and pi = -dH/dSi or, in view of (|?7jl . 



dS,. 



(28) 



This equation can be extended to take dissipation into 
account by adding the friction force (j^SJ: 



Pi 



m 

' dS, 



(29) 



In this paper we mainly consider a system of two iden- 
tical qubits, then the subscript i can be dropped in the 
parameter notation: 



Ci — C, Li L, /oi — /q, Ri — R 



(30) 



(the external flux 0j, which is the control parameter in 
experiments, is assumed to be generally different for the 
two qubits). Then Eqs. and (|2ZI) become 



H 



pi 



p\ 



Kp\P2 



2(1-)- 

P% + CPj 
(1 + Oim 



Ui{Si) + U2iS2 



(31) 
(32) 



where C = Cx/{C+Cx)- Experimental parameters of Ref. 

= 6 fF and Tl = 25 ns, correspond to C = 8.5 x 10"^ 
and R =35.7 kfl. 

The above formalism provides a basis for the analysis 
of the qubit-system evolution before the measurement 
(the gate operation) , during the measurement pulse, and 
after the measurement pulse (the crosstalk). The latter 
stage is considered below. 
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III. MEASUREMENT CROSSTALK: FIRST 
QUBIT DYNAMICS 

In the fast measurement scheme employed in Refs. 'g' 
and0, a short flux pulse applied to the measured qubit 
decreases the barrier between the two wells (see Fig. ^ , 
so that the upper qubit level becomes close to the barrier 
top. In the case when level |1) is populated, there is a 
fast population transfer (tunneling) from the left well to 
the right well. Due to dissipation, the energy in the right 
well gradually decreases, until it reaches the bottom of 
the right well. In contrast, if the qubit is in state |0) 
the tunneling essentially docs not occur. The qubit state 
in one of the two potential minima (separated by almost 
<l>o) is subsequently distinguished by a nearby SQUID, 
which completes the measurement process. 

In a system of two identical coupled qubits, crosstalk 
can produce measurement error if the qubits are in dif- 
ferent logical statesi^ For definiteness, we assume that 
before the measurement the qubit system is in the state 
1 10), i.e., the first qubit is in the excited state and sec- 
ond qubit is in the ground state. Then after the mea- 
surement the first qubit performs damped oscillations in 
the right well, which in the classical language^'' produces 
an oscillating (microwave) voltage {^q/2tt) 6i{t). This 
voltage causes oscillating current through the coupling 
capacitor Cx, which perturbs the second qubit. The ef- 
fect is so strong that measurement of the second qubit 
after the dissipative evolution of the first qubit is practi- 
cally useless: there is a little chance for the second qubit 
to remain in the ground statei^ The effect of crosstalk 
can be significantly suppressed if the two qubits are 
measured practically simultaneously (experimentally, not 
more than few nanoseconds apart) because the crosstalk 
excitation of the second qubit takes finite time>2i Never- 
theless, crosstalk leads to noticeable measurement errors 
even in the case of simultaneous qubit measurement. The 
reason is that strong excitation of the second qubit may 
lead to its switching from the left to the right well even 
in absence of the measurement pulse. 

This is exactly the effect which we analyze in this pa- 
per. We assume that the first qubit is switched to the 
right well at t = 0, while the second qubit at this time is 
in the ground state and no measurement pulse is applied 
to the second qubit (physically, this means that the pulse 
is short and does not change the qubit state). Our main 
goal will be analysis of the measurement error, which in 
this case is switching of the second qubit to the right well 
due to the crosstalk excitation. 

A rigorous theoretical study of the measurement 
crosstalk should involve a numerical solution of quan- 
tum evolution of two coupled qubits with the account 
of dissipation, which would require extensive computer 
resources. In the present paper we employ several sim- 
plified approaches, which have the advantage of being 
relatively fast numerically, thus facilitating a study of 
the crosstalk dependence on the parameters. The first 
qubit is always treated classically, while the second qubit 



is studied both classically and quantum mechanically. 

We will mainly consider two experimentally relevant 
cases of the second qubit biasing, characterized by the 
dimensionless barrier heights N12 = 5 or 10 (02 = 5.09 
and 4.84, respectively). Then the plasma frequencies are 
u;i2/2tt = 8.91 GHz and 10.2 GHz, respectively. The 
crosstalk mechanism is obviously very efficient when the 
first qubit oscillation is in resonance with uji2. 

For the first qubit we choose the biasing parameter 
Nil = 1.355 (01 = 5.31) which is close to the experimen- 
tal value^-^^ at which the state |1) efficiently tunnels out. 
The corresponding WKB tunneling ratc^'^ of the state 
1 1) is 3 X 10^ s~\ which ensures tunneling during few- 
nanosecond-long measurement pulse as in experiments 
of Refs. .6 and g. Notice that the barrier height Nn is 
smaller than the naive estimate 1.5 for the dimensionless 
energy of the state |1). Actually, because of significant 
anharmonicity, the energies of states |0) and |1) in this 
case are 0.475 and 1.26 from the well bottom in units 
of Ticuii, where uiii/2tt ~ 6.87 GHz. We neglect the fact 
that in the experiment, after the measurement pulse, the 
biasing of the left qubit returns back to Nn ~ 5, because 
as we checked, this does not lead to a significant change 
of the evolution dynamics in the right well. 

At the initial moment t = Q the first qubit is assumed 
to be in the right well close to the barrier top position, 
with the velocity 5i = 0. However, instead of assuming 
its initial energy to be the same as the energy of state 
|1), we choose a slightly lower energy, which is below the 
top of the barrier by 20% of the well depth A?7;i. This 
(somewhat arbitrary) choice prevents unphysically slow 
dynamics of a classical particle in the case when it is 
very close to the barrier top (quantum dynamics due to 
dissipation does not significantly slow down at the energy 
close to the barrier top). We have also checked that the 
qubit dynamics is not too sensitive to the choice of initial 
energy (when it is above the left well bottom and not too 
close to the barrier top). 

Since the initial energy of the first qubit with respect 
to the bottom of the right well is much higher than the 
maximal energy of the second qubit in the left well, one 
can neglect the back action of the second qubit onto the 
first one while the the second qubit remains in the left 
well. Then the second term in the numerator in Eq. H32|) 
for the first qubit (i = 1) can be dropped, yielding 



5i = pi/m" , 



(33) 



so Eq. H29|l with i — 1 gives the approximate equation of 
motion for the first qubit: 



5i 



5i 



27r/o 



C"R $oC"' 



sin Si 



Si - 01 
C"L 



0, 



(34) 



where C" = (1 + C)C. This is obviously the usual equa- 
tion for an isolated first qubit with capacitance C re- 
placed by effective capacitance C" = C + 0^0/ {C^ + C) 
which takes into account the series connection of the cou- 
pling capacitance and the second junction capacitance 
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(this corresponds to the approximation of zero charge at 
the second qubit). 

Note that even though the set of equations of motion 
(|18|l is equivalent to Eqs. H26|) and 129|) . the above ap- 
proximation makes them different. In particular, in the 
case of identical qubits the equation for the first qubit 
obtained from Eq. H18(l by neglecting the r.h.s., differs 
from Eq. H34|) by the substitution C" C . However, 
for small dimensionless coupling C, (which is the experi- 
mentally relevant case assumed here), the two equations 
differ by very small terms on the order of C^- Physically, 
C = C + Cx as the effective capacitance of the first qubit 
corresponds to the model in which the voltage across the 
second junction is neglected (in contrast to the charge in 
the previous model). 

Equation (|34|l shows that the first qubit performs 
damped non-harmonic oscillations. Because of anhar- 
monicity, the gradual decrease of the qubit energy Ei = 
rn! -|- C/i(^i) due to dissipation leads to the grad- 
ual increase^ of the oscillation frequency (driving the 
second qubit) which can be obtained asii 



f^\E,) = V2^ / , (35) 



where a and b are the classical turning points. The time 
dependence fd{t) of the oscillation frequency is shown in 
Fig. for Cx — (solid line) and 6 fF (dashed line) as- 
suming Ti — 25 ns. The curves are very close to each 
other showing very small effect of the capacitance renor- 
malization (so that even smaller second-order effect of 
the difference between C" and C" is really negligible). 
We have also checked numerically that variation of Ti 
from 25 ns to 500 ns does not change noticeably the de- 
pendence fd{t) if the time is normalized by Ti, which is 
rather obvious since fd ^ Tf^. 

Figure |31 shows that the oscillation frequency sharply 
increases initially and then slowly tends to the right-well 
plasma frequency ujri/^ir ~ 15.3 GHz (the dash-dotted 
horizontal line in Fig. O . This is explained by the fact 
that the initial system energy is close to the barrier top, 
where the oscillation frequency is significantly lower (it 
tends to zero when the energy approaches the barrier 
top), while anharmonicity becomes relatively weak after 
the energy is no longer close to the barrier top. 

Notice that it takes a finite time for the first qubit 
dynamics to get into resonance with the second qubit 
((JJ12I2-K is around 9-10 GHz, as mentioned above); we 
find from Fig. 13 that = 0.085 Ti for = 5 and = 
0.192 Ti for N12 ~ 10. As a simple estimate, this is the 
time after which the second qubit becomes significantly 
excitediSi 




0.2 0.4 0.6 0.8 1 1.2 1.4 
tITi 

FIG. 5: The first-qubit oscillation frequency fd as a function 
of time t (normalized by the energy relaxation time Ti) for 
Ci^ = (solid line) and Ca; = 6 fF (dashed line), assuming 
Nil = 1.355 and parameters of Eq. ^V2\. Dash-dotted hor- 
izontal line, Wt-i/Stt = 15.3 GHz, shows the long-time limit 
of fd{t). Two dotted horizontal lines show the plasma fre- 
quency for the second qubit: u}i2/2n = 10.2 GHz for N12 = 10 
and ijji2/2n — 8.91 GHz for N12 = 5. The arrow shows the 
moment tc of exact resonance in the case N12 = 5. 



IV. SECOND QUBIT: CLASSICAL APPROACH 

In the approximation (|33() . the second-qubit equation 
of motion follows from Eqs. H32(l and 1)29(1 with i — 2: 

^.+ (Y^ + ^Bin^. + -^ = C^,(.). (36) 

Here we assume that the second-qubit relaxation time 
T[ may be different from the first-qubit value Ti. For 
simplicity we will mostly neglect the energy relaxation 
for the second qubit (r{ = 00), while more experimen- 
tally relevant case T{ — Ti will be discussed only in the 
classical approach. 

Equation l|36|) has a simple physical meaning as an 
evolution of the second qubit with effective junction ca- 
pacitance C", externally driven by the oscillating current 
{^o/2tt)(C"5i. However, considering oscillating voltage 
(<i>o/27r)5i across the first junction coupled to the second 
qubit via capacitance Cx, one would expect the driving 
current to be ($o/27r)Cx^i [as in Eq. lfTF|l ]. The relative 
difference between ^C" and Cx is on the order of <C 1, 
which is the accuracy of treating two qubits separately. 

In this section we consider classical dynamics of the 
second qubit due to the oscillating driving force iii 
Eq. H3t)f) with a slowly varying period and amplitude. 



A. Harmonic-oscillator model 

Though Eq. H36|l is simpler than the exact equations 
(|22Jl and (|25J| , the behavior described by Eq. is still 
complicated and generally chaotic. To get an insight 
into the behavior of the second qubit, we first model it 
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by a harmonic oscillator )S i.e., substitute Eq. H36|l by 

x + oj^^x = Ch{t), (37) 

where x = 82 — 5i2, S12 is the left-well minimum position, 
and is the unperturbed plasma frequency. Actually, 
the small-vibration frequency in Eq. should be dif- 
ferent from UJ12 due to the mass (capacitance) renormal- 
ization m m" , so that CO12 uJn/y/l + C', however, 
for small coupling considered here {( < 1%) we neglect 
the difference. Correspondingly, the oscillator energy is 
E2 = m{x^ +10^2^"^) 1"^ (here we also neglect the difference 
between m and m"). 

We have simulated Eq. H37() numerically, assuming the 
system initially at rest at the potential minimum, a;(0) — 
i(0) = 0, and checked that the oscillator energy coincides 
with the analytical solutionsii 



E2{t) 



t 2 
dt'e-'"'^'"5i{t') . (38) 



Figure shows the time dependence of the energy E2 (t) 
in units of hLUi2 for uji2/2Tr — 8.91 GHz, which corre- 
sponds to N12 — 5 (parameters of the first qubit evolu- 
tion have been discussed in Sec. Ill and correspond to 
Fig. |SJ). One can see that the energy E2 remains very 
low until a sharp increase followed by gradually decreas- 
ing oscillations. This behavior can be easily explained by 
changing in time frequency /d of the driving force (Fig. 
O which passes through the resonance with the second 
qubit (2i 

For an analytical analysis, let us consider the vicinity 
of the moment of exact resonance, Wd(ic) — ^12 (here 
Lo^ = 271 J d), and approximate 5i{t) as a harmonic signal 
5i{t) = A{t) exp(i J ujd dt) with constant amplitude A = 
A{tc) and linearly varying frequency 



UJd{t) = iOl2 + ct{t - tc), 



(39) 



with a — ujditc) > (we have neglected the complex- 
conjugated term and higher-order harmonics as being out 
of resonance). Thus the problem is reduced to passage 
of a harmonic oscillator through resonance with a con- 
stant ratefSfliSi Also assuming slow crossing, a <C ujf2, 
and shifting the lower endpoint of integration in Eq. 1381) 
to —00 (which is a good approximation for tc 3> y/2/a), 
we obtain 



E2{t) = EoF{i), Eo 



F{i) = 



t 



t-tr 



/ e'^^'di] 


2 

1 




^ 4 




J —00 







(40) 



Notice that the function F{t) with i proportional to a 
spatial coordinate describes the Fresnel diffractionSi and 
has the following asymptotic dependence: 



F{i) 
Fit) 



1 



sin(P - 7r/4) 



/TTt 

for 



fori>l, (41a) 
i>l. (41b) 
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FIG. 6: The second qubit energy E2 (in units of huji2) in 
the oscillator model as a function of time t (in ns) for (a) 
C:, = 5 fF and Ti = 25 ns and (b) = 2.5 fF and 5 £F 
and Ti = 500 ns, while N12 — 5. Dashed line in (a) shows 
approximation using Eq. H4()^ . The arrows show the moment 
tc when the driving frequency fa (see Fig. |5J| is in resonance 
with uji2/2-n- = 8.91 GHz. 



The oscillating term in Eq. (|41a|l describes the beating 
between the oscillator and driving force frequencies, with 
the difference frequency increasing in time, d(t'^)/dt = 
ot[t — ic), and amplitude of beating decreasing as 1/i 
(see dashed line in Fig. Et). Notice that i^(0) = 1/4, 
F{po) = 1, and the maximum value is F(1.53) = 1.370, 
so that Eq is the long-time limit of the oscillator energy 
i?2, while the maximum energy is 1.37 times larger: 



Eo 



\.?,lT:CmA^ 



(42) 



This result for the energy can be understood in the 
following way. In case of exact resonance, the oscillation 
amplitude of x(t) in Eq. H37|) increases linearly in time 
with the rateii C^AIijJi2. The effective time of resonance 
At corresponds to a significant phase shift due to beating: 
Q!(A<:)^ ~ 1. Therefore the resulting amplitude after the 
resonance crossing is '-^ C,A/uji2\/a and the corresponding 
energy is E2 ^ (^mA^ ja^ which coincides with Eq. \\2\ 
up to a numerical factor. 

For the parameters of Fig.|Sl(iV;2 — 5) we find from Fig. 
ISjthat exact resonance between /d(t) and w/2/27r = 8.91 
GHz occurs at t^ = 0.085 Ti (in particular, tc = 2.13 ns 
for Ti = 25 ns and t^. = 43 ns for Ti = 500 ns) and 
a = 110ns~VTi. To compare Eq. (gUIl with the nu- 
merical results, we also need the value of A. It can be 
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estimated a.s A = ujf^^/'^^ where A is the ampHtude of 
Si{t) oscihations at t — (the factor 1/2 comes from 
our definition of A as half of the amphtude of 5i oscilla- 
tions). Using the numerical result A ~ 2.7 corresponding 
to the resonance with the second qubit at N12 = 5, we 
find A = 4.3 x lO'' ns^^ (this resonance happens at 121 
GHz below the barrier top; A is defined as half of the full 
span oi Si oscillations - see Fig. The dashed line in 
Fig. Ela) shows the corresponding analytical result (|42|l . 
One can see that the analytics fits the oscillator energy at 
t = tc pretty well; however, the maximum energy -E2,max 
given by Eq. (|42|l is somewhat different from the numer- 
ical result: -B2,inax/?iw;2 = 11.8 versus 14.2 numerically. 
(The difference decreases with increase of Ti : for Cx = 5 
fF and Ti = 500 ns the corresponding numbers are 235 
and 242.) A noticeable discrepancy between the analyt- 
ical and numerical results in Fig. Efa) can be attributed 
mainly to the fact that A and a change with time, in 
contrast to the assumptions made in the derivation of 
Eq. H4()|l. It is interesting to notice that A{t) initially 
increases because of the frequency fd increase, while it 
starts to decrease at i > 0.52 Ti (after reaching the max- 
imum of 5.8 X 10"^ ns~^) because of A decrease. Actually, 
a good fit of the numerical results by the values of A 
obtained as ^ = (27r/d)^A/2 is quite surprising, because 
oscillations Si (t) are strongly non-harmonic (even having 
three maxima and three minima per period). We have 
also calculated A using some advanced analysis of Si{t) 
and found values very close to the simple estimate. (This 
other method is based on calculating Fourier transform 
of Si within a time interval around tc, cutting off the 
spectrum above the minimum between the first and sec- 
ond Fourier peaks, calculating inverse Fourier transform, 
and finding the oscillation amplitude at t = tc.) 

Increase of N12 leads to more efficient excitation of the 
second qubit. For example, for N12 — 10 and other pa- 
rameters as in Fig. |Hfa), the numerical maximum value 
of E2,niax/^0Ji2 bccomcs 30.5 (more than twice larger 
compared to the case = 5). This happens because 
of the decrease of a{t) and increase of A{t) with time, 
and correspondingly with Ni2. [For N12 — 10 (so that 
uji2/2n ^ 10.2 GHz), we find = 0.192 Ti (i.e., = 4.8 
ns for Ti = 25 ns and = 96 ns for Ti = 500 ns), 
a = 57ns-V7'i. and A{tc) = 5.2 x 10^ ns^^.] 

So far we have completely neglected the energy re- 
laxation in the second qubit. Since the effective time 
of resonance is At ~ S/^/a [rise time of the function 
F{i) from the 10% level to the maximum], the neglected 
effect should not be important (less than ^^10%) for 
Ti ^ 30/y/a. Using the estimate a ~ llOns^^/Ti (see 
above), we find that taking into account the second qubit 
relaxation would not change significantly our results for 
-E'2,max ff Ti JSJ 10 ns, which justifies our model. 

Let us discuss the dependence of the maximum energy 
-E'2,max of the second qubit on C'x and Ti. Taking into 
account that ( (X Cx (for Cx ^ C) and a oc 1/Ti, we 



obtain from Eq. H42|l the scaling 

S2,max (X C^Ti. (43) 

As seen from Fig. El numerical results confirm the obvious 
scaling i?2,max C!C Cx, while the scaling i?2,inax c>c Ti is not 
very accurate, but is still good as a first approximation. 

In this subsection we have treated the second qubit 
as a harmonic oscillator. However, to analyze the mea- 
surement error due to the crosstalk, we have to as- 
sume switching from the left well to the right well when 
£-2, max > Ni2hu;i2 (which is surely not fully consistent 
with the oscillator model). All curves in Fig.|Hl(A'/2 — 5) 
correspond to such switching, leading to the measure- 
ment error. The measurement error can be improved by 
decreasing the coupling capacitance Cx , which should be 
chosen to be smaller for larger Ti . Equation H43|l implies 
that to avoid the errors due to crosstalk, one needs to 
choose 

Cx < Cx^T = B/^i, (44) 

where Cx,t is the threshold coupling capacitance. From 
the numerical simulations (see Fig. EJ we obtain S ~ 15 
fFy^ in the case Na — 5. For N12 = 10 we get 
_B ~ 14 fF -yns, yielding a little stricter condition than 
for Nn = 5. Notice that for experimental parameters of 
Ref.|9r(C:r=6 fF, Ti =25 fF) this bound is exceeded ap- 
proximately twice, which is an indication that our simple 
model is not sufficiently accurate. As we will see in the 
next subsection, the theoretical bound is softer (higher) 
when we use actual potential profile for the second qubit 
instead of using the harmonic oscillator model. 

B. Actual qubit potential 

Let us analyze the second qubit evolution still using a 
classical model, but taking into account the actual po- 
tential profile 1/2(^2), i-e., solving Eq. instead of the 
simplified equation (|S7|l . Figures [7| and |H1 show the time 
dependence of the second-qubit energy, 

E2^m"52^ /2 + U2{S2), (45) 

in the absence of dissipation in the second qubit {T[ — 
00) for A^;2 = 5 and 10, while Ti = 25 ns. (In this sub- 
section we take into account the mass renormalization 
m — > to" explicitly, even though this does not lead to 
a noticeable change of results.) A comparison of Figs. 
EJa) and [3 shows that in both models the qubit energy 
remains small before a sharp increase in energy. How- 
ever, there are significant differences due to account of 
anharmonicity: (a) The sharp energy increase occurs ear- 
lier than in the oscillator model (the position of short- 
time energy maximum is shifted approximately from 3 
ns to 2 ns); (b) The excitation of the qubit may be to a 
much lower energy than for the oscillator; (c) After the 
sharp increase, the energy occasionally undergoes notice- 
able upward (as well as downward) jumps, which may 
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FIG. 7: The second-qubit energy E2 (in units of hun) as a 
function of time t (in ns) for N12 — 5, Ti — 25 ns, and (a) 
= 5 fF, (b) = 6 fP, and (c) C^, = 9 fF. The classical 
model with actual qubit potential is used; energy relaxation 
in the second qubit is neglected, Ti = 00. In (b) and (c) the 
qubit switches (goes over the barrier) at 44 ns and 2.1 ns, 
respectively. 



overshoot the initial energy maximum; (d) The model 
now explicitly describes the qubit escape (switching) to 
the right well [Figs.[7Jb) and[7{c)]; in contrast to the os- 
cillator model, the escape may happen much later than 
initial energy increase; for example, in Fig. [TJ^b) the es- 
cape happens at t ~ 44 ns ^ — 2.1 ns. 

The properties (a) and (b) can be understood by taking 
into account the fact that the oscillation frequency in the 
second qubit decreases with the energy increase (it should 
become formally zero at the top of the barrier), while the 
driving frequency increases with time (Fig.lsj. Therefore, 
initially small out-of-resonance beatings when ujd < LOn 
are amplified because of the positive feedback: larger am- 
plitude makes it closer to the resonance, which increases 
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FIG. 8: Same as in Fig. □ for Nn = 10, Ti = 25 ns, and 
C:. = 6 fF. 



the amplitude even more. This makes the non-excited 
state unstable, which leads to a sharp increase of the 
qubit energy earlier then the condition LOd = ^12 is sat- 
isfied. The same mechanism is also responsible for lower 
qubit excitation, when compared to the harmonic oscil- 
lator model: the resonance cannot be as efficient as in 
the harmonic oscillator model since the qubit excitation 
quickly moves the qubit frequency out of the resonance. 
The property (c) is related to crossing of higher-order 
resonances, which occur when tUdit) is commensurateii 
with the oscillation frequency of the system, which itself 
depends on the energy E2{t) and hence on the time. Sim- 
ilar mechanism is responsible for the qubit switching at 
t ^ tc', in particular, in Fig.[7Jb) the switching happens 
when the driving frequency fd becomes approximately 
twice larger than the second qubit frequency. 

In contrast to the oscillator model, Eq. (|36() for the ac- 
tual qubit potential cannot be solved analytically,?- so we 
rely only on the numerical simulations. We are interested 
in the conditions, for which the system remains in the left 
well. Generally, the qubit excitation increases with in- 
crease of the coupling Cx ; therefore one expects a certain 
threshold value Cx.t (depending on Ti and other param- 
eters), which separates the switching and no-switching 
scenarios. However, because of the complex dynamics of 
the system, the dependence on Cx is non-monotonous, 
so that increasing Cx may sometimes change switching 
case into no-switching case. In this situation, we define 
Cx,T as a minimum value at which the switching happens 
(even though larger Cx may correspond to no-switching) . 
Similar to the harmonic oscillator model, we expect that 
Cx.t generally decreases with increase of Ti ; however, be- 
cause of the complex dynamics, the dependence Cx,t(Ti) 
should not necessarily be monotonous. 

The dots connected by two solid lines in Fig. |51 show 
the numerically calculated Cx.t for N12 — 5 and 10, for 
5 values of Ti ranging from 25 ns to 500 ns (so far we 
still assume T{ — 00). For these calculations we have 
used the increment of 0.1 fF for Cx and simulated the 
qubit dynamics in the time interval < i < 6T1. Notice 
that the lines are not smooth (the lower line even has a 
bump), which is the result of irregular nonlinear dynam- 
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FIG. 9: The log- log plot of the threshold coupling capacitance 
Cx,T vs. Ti in the classical model assuming = oo (solid 
lines) or T{ = Ti (dashed line) for N12 = 5 (dashed and lower 
solid line) and 10 (upper solid line). The crosstalk excitation 
does not switch the second qubit if Cx < Cx,t- The numerical 
data are shown by the dots; the lines are just guides for the 
eye. 



ics of the system. Nevertheless, the numerical results 
confirm the generally decreasing dependence Cx,t(Ti), 
and fitting solid lines by the formula 



(46) 



(where Cx,t is measured in fF while Ti is measured in ns) 
we obtain i? ~ 8 and 12 for N12 = 5 and 10, respectively, 
while (3 ~ 0.12 for the both lines. 

Since the gate speed is proportional to Cx, it is ad- 
vantageous to have higher Cx- The above results show 
that raising the barrier after the measurement pulse to 
N12 = 10 would allow us to increase Cx i^ comparison 
with the case N12 = 5. The reason for this can be under- 
stood by comparing Figs. [7{b) and|Sl which show that 
in the case N12 — 10 the sharp energy increase is lower 
relative to the barrier top than for N12 = 5. Note that 
the dependence (|46|l is much weaker than the relation 

— 1/2 

Cx.T oc T-^ obtained in the oscillator model, which is 
advantageous for design of qubits with weak decoherence 
(large Ti). 

Now let us consider the effect of dissipation in the sec- 
ond qubit. The dashed line in Fig. O shows Cx.t{Ti) 
dependence in the presence of dissipation (r{ — Ti) for 
the same other parameters as for the lower solid line 
(for which T{ = 00). As we see, account of dissipa- 
tion increases Cx.t quite noticeably, which contradicts 
the conclusion from the harmonic oscillator model (pre- 
dicting no significant dependence) . The reason is that for 
Cx slightly above Cx,t the switching in the model with- 
out dissipation usually occurs significantly later than the 
initial sharp increase of the energy [see Fig. [7[b)] and 
is caused by "secondary" jumps of the energy due to 
strongly nonlinear dynamics, as discussed above. Dis- 
sipation in the second qubit (Fig. I10|) shortens signifi- 
cantly the time interval during which the switching due 
to secondary jumps is possible, thus increasing Cx,t- (We 




FIG. 10; The second-qubit energy E2{t) in the classical model 
taking into account energy dissipation in the second qubit, for 
N12 = 5, = 6 fF, and Tl = Ti ^ 25 ns. [Compare with 
Fig.Hb).] 



would also like to mention a possibility of a system return 
into the left well after the escape into the right well, which 
may take place with or without dissipation.) Fitting the 
dashed line in Fig. 1^1 by the power-law dependence 
we find B ~ 12 and (3 ~ 0.13, so that the scaling power 
(3 is practically the same as in the no-dissipation case, 
while B becomes considerably larger. 



V. SECOND QUBIT: QUANTUM APPROACH 

In the quantum approach the second qubit is described 
by the wave function ^'((5, t) (in this section we often omit 
the subscript 2 to shorten notations), which obeys the 
Schrodinger equation 



Here the Hamiltonian 



H{t) 



2(1 + C)m 



UiS), 



(47) 



(48) 



in which p = —ifi{d/dS), follows from Eq. H31|l by con- 
sidering p2 as the operator p, while pi (t) and Si (t) are 
considered as classical functions of time obtained from 
Eqs. (|34|l and H33|l : the first-qubit energy in this case 
does not contribute to the Hamiltonian (|48|l . 

The term linear in p in Eq. (|48|l has the same form 
as for the interaction of a charged particle with a time- 
dependent electric field described by a vector potentialiS^ 
Using the gauge transformations^ 



*((5,t) = *'((5,i)e-^^fi(*)'^/'*. 



(49) 



we can replace the vector-potential by a scalar potential 
in the Hamiltonian. Then Eq. H47|l becomes 



ih 



'dt 



H'{t)^', 



(50) 
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in which the Hamiltonian [subtracting c-number term 
C^pl{t)/{2m") and using Eq. becomes 

H'it)=Ho + Vit), 

Ho = ^ + U{S), V{t)=-Cm%it)S. (51) 

This Hamiltonian exactly corresponds to the classical 
model used in the previous section. Similarly to the 
classical case, the difference between Cm" and rrix = 
{^o/2tt)'^Cx in the formula for V{t) should not be taken 
seriously (as being within the accuracy of treating two 
qubits separately). 

The partial differential equation l|5U|l can be reduced to 
an infinite set of ordinary differential equations jSS. using 
the expansion of the wavefunction over the eigenfunctions 
MS) of Ho, 

*'((5,i)=V Unit) MS)- (52) 

^ — 

Inserting expansion 1)52(1 into Eq. H50|l yields the set of 
equations for the coefficients a„(t): 

a,i = -l(A„/n)a„ H > Onn'On', {o3) 

where i?„ is an eigenvalue of and bnn' is the "position" 
matrix element: 

/C30 
MS)SM{S)d5 (54) 
-oo 

(notice that we will use a different notation for the Kro- 
necker symbol). 

We have calculated the eigenstates and eigenvalues 
of Hq numerically, using the Fourier grid Hamiltonian 
method^^ (same as periodic pseudospectral method^^). 
After obtaining eigenfunctions, we have calculated the 
matrix (5„„' (|54|l and solved numerically Eqs. (|53|l . re- 
stricting the space to a finite subset of the states. For 
a given 7V;2, we take a reasonably small number of 
consecutive states (n = n^, n.^ + 1, . ..,71^ + 71^), which in- 
clude all left-well states and provide a sufficiently good 
approximation to the exact solution (the choice of the 
subset of states is discussed below). The column vec- 
tor a — (a„;, . . . ,a„;+„^)^ satisfies the equation [cf. Eq. 

m] 

ind = H{t)a, (55) 

in which the matrix 

Hnn'{t) = EnS^^, - Cm"Si{t)Sn„' (jT-i <n,n' < m + rir), 

(56) 

is the Hamiltonian of the system in the restricted Hilbert 
space spanned by the subset of states (here (5^, is the 
Kronecker symbol). 

We define the probability Pi (t) to find the second qubit 
in the left well as 

m)=J2 Prnit). Pn{t)^\an{t)\\ (57) 
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FIG. 11: The second qubit potential U{5) (thick line) and 
eigenfunctions |'!/'n(5)p (145 < n < 171), shifted vertically by 
the energy eigenvalues En (in units of Tiuan) for Nn = 10 and 
Ci^ = 6 fF (tj;2/27r = 10.2 GHz). The energy origin is chosen 
at the minimum of the left well. An arbitrary scale for \ipn\'^ is 
chosen to be 1/8 of the E„/hLJi2 scale. The arrow illustrates 
the dominating transition responsible for the escape from the 
level fc = 9 at t ~ 16 ns. 



where Pn{t) is the probability of state n occupation, and 
the summation is only over the states localized in the 
left well. We will also denote P„, (t) as Qk{t), where 
k enumerates the states in the left well, starting from 
k = (the left-well ground state). We define the switch- 
ing probability as Ps{t) = 1 — Pi{t). Notice that we 
consider transition to delocalized states (above the bar- 
rier) as escape from the left well (even though in this case 
there is a possibility of "repopulation" of the left well if 
dissipation is taken into account). 

Figure shows (for the case N12 = 10) the eigen- 
functions |'0„((5)p and the corresponding energies En for 
145 < n < 171, where the level numbering starts with 
n = for the ground state (in the right well). One can 
distinguish 3 types of states: (a) 12 states localized in 
the left well (n = 146, 148, . . . , 166, 169 or, respectively, 
k = 0,. . . , 11), (b) delocalized states (n > 170), and (c) 
states localized in the right well (the remaining states). 
Depending on the barrier height N12 (controlled by the 
external flux 02), the resonant states may also be present: 
when the energies of states localized in the left and right 
wells approach each other sufficiently close, the states 
mix and become delocalized. Actually, in Fig.^Jthe left- 
well states A: = 10 and 11 are partially delocalized due to 
interaction (tunneling) with neigboring right-well states. 
Notice that even though the left-well energies are practi- 
cally insensitive to the coupling capacitance Cx {C'x — 6 
fF for Fig. nil) , their relative energy shift with the right- 
well level comb depends on Cx significantly because the 
right well is very deep, n > 10^. 

FigurelT^shows the left- well population Pi (t) for A'';2 — 
10 and Ti — 25 ns (energy dissipation in the second qubit 
is neglected, T{ = 00). The subset of the states used 
in this calculation is 141 < n < 185. Figure [TSI shows 
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FIG. 12: The probability P; = 1 — Ps for the second qubit to 
remain in the left well as a function of time t for N12 — 10, 
= 6 fF, and Ti = 25 ns. 



the populations of the first 10 levels in the left well; the 
populations of levels fc = 10 and 11 are not shown since 
they are very close to zero. 

Figure lTH shows the time dependence of the qubit mean 
energy 



(58) 



for the same parameters as in Figs. El and El Com- 
paring the mean energy with the classical qubit energy 
for the same parameters (Fig. jSl, we see that the two 
curves are similar, though classical energy shows larger 
fluctuations. Note that the mean energy starts at i = 
from a nonzero value equal to the qubit energy in the 
ground state « huji2/2. Even though the mean energy is 
significantly lower than the barrier height (similar to the 
classical energy), the escape probability Ps{t) = 1 — P;(t) 
is nonzero in the quantum case (see Fig. lT^ . 

The time dependence of the switching probability (Fig. 
II2II looks quite irregular. This shows that the quantum 
behavior is rather complicated, similar to the classical 
behavior discussed in Sec. IIVBI As seen in Figs. 1131 and 
1141 at t < 3 ns, when the driving force is far from the 
resonance with the qubit, the population mainly remains 
in the ground state. At this stage, there is no switching 
(see Fig. ll2|) . Similar to the classical case, there is a fast 
qubit excitation (though still almost without switching) 
between 3 ns and 4 ns (a little earlier than the moment 
tc — 4.8 ns of classical resonance), while the main switch- 
ing happens much later, mostly at < ~ 16 ns. 

To understand the excitation mechanism, we show in 
Fig. Elthe Rabi frequencies Rk,k-i — 2CnT-"^|<5/c,/c-i 
for the adjacent left-well transitions; Rk,k-i is equal 
to the amplitude of Vk,k~i{t) /fi oscillations [actually, 
since V{t) oc Siit) given by Eq. (|51|l is significantly 
non-harmonic in time, we need to use the amplitude 
of the resonant component]. Since the amplitude 2A{t) 
of 5i oscillations changes with time, the Rabi frequen- 
cies also change with time. In Fig. El we show the 
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FIG. 13: The occupation probabilities Qk(t) for the left- well 
levels k as functions of time t (in ns) for A^; = 10, d = 6 fF, 
and Ti = 25 ns. 



values corresponding to the exact classical resonance, 
fditc) = i^vz/^TT, which happens at tc = 4.8 ns (then 
A — 5.2 X lO'^ ns^^); for comparison, at i ~ 3 ns the 
value of A and, correspondingly, the Rabi frequencies 
are approximately 10% smaller. (Notice that the mo- 
ment tc scales with Ti , but the values of Rabi frequen- 
cies at tc do not change with Ti.) For the levels not 
too close to the barrier top, one can use the harmonic- 
oscillator relation^^ |<5A;.fc-i| ~ \/kfi/2m"u}i2 ~ O.IOa/A: 
[for N12 = 10 and parameters of Eq. H12(l ]. yielding 
Rk,k-i/2n w A{C^/C) ^y2km"/huJi2/2TT ~ l.lVk GHz. 
This formula fits well the numerical results in Fig. ll5l up 
to fc = 9; for higher levels anharmonicity becomes really 
strong. 

Figure 1161 shows the time dependence of the corre- 
sponding detunings a;fe^fe_i/27r — fdit)- Though the ex- 
act resonance with the transition <-> 1 (cLiio/27r = 10.0 
GHz) happens att = 4.3 ns, a significant excitation starts 
earlier, at i ~ 3 ns, when the detuning cjio/27r— fdit) be- 
comes less than the Rabi frequency Riq ~ 1 GHz. Since 
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FIG. 14: Time dependence of the mean energy (E) of the 
second qubit [Eq. H58|l ] in units of hui2 for N12 = 10, Cx — 6 
fF, and Ti = 25 ns. 
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FIG. 16: Time dependence of the frequency detunings 
[u>k^k-i — i^d(i)l/27r of left-well transitions {k — 1) ^ k for 
A''; = 10, Cx = 6 fF, and Ti = 25 ns. Curves from top to 
bottom correspond toA: = l,2...11. 
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FIG. 15: Dots: Rabi frequencies Rk^k-if^ir for the left-well 
transitions at f = tc, for Ni = 10, Cx = 6 fF, and Ti = 25 ns. 
Dashed line shows analytical dependence 1.1%/^ GHz. 



the Rabi frequency increases with fc, while the detuning 
first decreases and then increases with k (after detuning 
changes sign) , the ground level population rapidly propa- 
gates to higher levels, until the detuning becomes so large 
that the further excitation stops. As a result, at t ^ A 
ns the levels 0-2 become practically empty, while almost 
all the population is transferred to levels 4-7. Similar to 
the classical case, the excitation efficiency is significantly 
suppressed by the fact that driving frequency fd{t) in- 
creases with time, while the level spacing decreases with 
the level number k. Therefore, by the time at which levels 
k }^ 5 become populated, the further up-transitions are 
already out of resonance (which happened for them ear- 
lier) , becoming even farther off resonance with increasing 
time. 

The escape (switching) in the quantum case can occur 
in several ways. The population which goes to the highest 
states A; = 10 and 11 is lost rather fast (within a fraction 
of nanosecond), since those states interact significantly 
with the right-well and delocalized states. However, in 
our simulation transitions to these states are quite weak 
because of significant detuning (even though the Rabi 
frequency is not much smaller than detuning, the level 
occupations above k = 7 decrease 4-10 times per level, so 



the occupation of the level /c = 10 is already very small). 

Another switching mechanism is the following. With 
the increase of the driving frequency fd{t), it can become 
resonant with transitions between non-adjacent states, 
thus populating the states close or above the barrier, 
which cannot be populated otherwise. In particular, 
more than half of the switching probability in Fig. 1121 is 
due to the sharp decrease of Pi {t) between 16 ns and 18 
ns, which happens because of the transition between the 
state k = 9 {n = 164) and the right-well state n = 168 
with the difference frequency of 13.2 GHz (see the arrow 
in Fig. Ill() . At i = 16 ns the detuning for this transition 
is 0.7 GHz {fd — 12.5 GHz), which is much smaller than 
other detunings between the levels of interest and is com- 
parable to the corresponding Rabi frequency ~ 0.5 GHz. 
This value of the Rabi frequency is obtained as above, 
taking into account the matrix element |(5i64,i68l = 0.048 
and the value ^ = 5.8 x 10^ ns~^ at t = 16 ns. 

Notice that in our numerical method the wavefunc- 
tion is represented in the basis of non-perturbed (time- 
independent) eigenstates. An alternative way would be 
to diagonalize the Hamiltonian at each moment of time 
and use the time-dependent eigenstates. Even though 
both methods are formally equivalent, the second method 
would be more natural to use if the dissipation is taken 
into account. In the time-dependent language an im- 
portant mechanism of escape is Landau-Zener tunneling 
through the barrier. The perturbation V{t) in the Hamil- 
tonian l|51|) is equivalent to changing in time magnetic 
flux. 



c^2-^c^2+CXm"Si{t)/Ej, 



(59) 



which changes N12 and leads to oscillations of the energy 
shift between the comb of levels in the left well and the 
right- well comb. Because of rather strong amplitude of 
these oscillations (in the example of Figs. I12H16I N12 os- 
cillates between 9.15 to 11.2), each left- well level crosses 
with several right-well levels during one oscillation cycle. 
These (avoided) crossings lead to transitions between the 
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states in different wells (tunneling), the rate of which, ac- 
cording to the Landau-Zener formular^ is oc , where 
W is the minimal level splitting at the crossing. The 
values of W increase exponentially with increase of k. 
Therefore, the Landau-Zener tunneling is a relatively 
slow (ineffective) switching mechanism for all levels, ex- 
cept for the highest ones (such as fc = 10 and 11). Besides 
the Landau-Zener mechanism, the escape from the left 
well may also happen because the upper left-well states 
may become delocalized (above the barrier) when the 
barrier N12 decreases in the process of oscillations. We 
would also like to mention that the oscillations of Si (and 
therefore of N12) are strongly non-harmonic. In partic- 
ular, the mentioned above range 9.15 < N12 < 11.2 re- 
mains constant during long time interval < i < 0.4 Ti, 
because these extrema actually do not correspond to the 
turning points of 6i(t) oscillations; instead, they corre- 
spond to the points of inflection S'^ and Sc + 27t of U (x) 
[see Eq. {TJ], as long as these points are within the oscil- 
lation swing of Si{t). 

We have performed extensive calculations of the 
switching probability Ps(t), varying the parameters Cx, 
Ti, and N12. We run simulations within the time interval 
[0,6Ti]; after 6T1 the first-qubit oscillations decay to a 
very low level, so the perturbation of the second qubit is 
weak, and the further change of Pg {t) is practically neg- 
ligible. Correspondingly, we define the total switching 
probability Pg (which is the crosstalk error probability) 
asP,(6ri). 

Before discussing the results, let us briefly discuss our 
choice of the restricted subset of states used for the 
numerical solution of Eq. (|55|) . First, we consider the 
eigenvalues of -ff„„' (|56|l as functions of Si, and require 
that this dependence is sufSciently close to the depen- 
dence obtained using the full Hamiltonian H51|l . Besides 
comparing the energies, we also compare the matrix el- 
ements Snn' obtained using either the full or restricted 
space (now only for the maximum and minimum val- 
ues of 5i), and require that the difference is below 1%. 
The discussed below calculations for the case N12 = 5 
have been performed for the subset of rir = 30 states: 
167 < n < 196 (in this case there are 6 states lo- 
calized in the left weU: n = 169, 171, . . . , 179). For 
N12 = 10, as mentioned above, we have used Ur = 45 
states: 141 < n < 185. We have also performed calcula- 
tions with enlarged subsets (up to rir = 150) and found 
that the switching probability varies irregularly with the 
number of states in the subset; however, this variation 
is not significant, at least in the case of low switching 
probability Pg ^ 0.3, which is the range of our interest. 

The results of numerical calculation of Ps{Cx) depen- 
dence for several values of Ti and two values of N12 (5 and 
10) are shown in Fig. El Notice that the lines connect- 
ing the data points are not smooth and sometimes are 
even nonmonotonous. This may be explained by compli- 
cated dynamics, similar to the irregular behavior in the 
classical case. Despite the Ps{Cx) dependence in Fig. 1171 
is not smooth, we still see that the switching probability 
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FIG. 17: The qubit switching (error) probability Ps as a func- 
tion of coupling capacitance Cx (in fF) for T\ — 25, 50, 100, 
200, and 500 ns for (a) N12 = 5 and (b) N12 = 10. The numer- 
ical data are represented by points, connected by solid lines 
as guides for the eye. The dashed straight lines are results of 
the least-squares fit (notice the logarithmic scale). 



Ps decreases approximately exponentially with decrease 
ofCx- 

Using the linear (on the semilog scale) least-square fit 
for the results shown in Fig. El (see dashed lines), we 
obtain the contour plots for Ps on the plane of Cx and 
Ti (see solid lines in Fig. [TH|l . The data in Fig. ITS] can 
be fitted by straight lines, yielding approximate power- 
law dependence [(similar to Eq. (|46|l ] for the threshold 
coupling capacitance Cx^t'- 



CxATi) ~ B{Ps)Ti 



(60) 



which now depends on the tolerable level Ps of the mea- 
surement error probability. The obtained numerical val- 
ues of the parameters B and /3 in this formula are shown 
in TableQ] Notice that (3 depends on Ps quite weakly, but 
decreases appreciably when Ni changes from 5 to 10. The 
values of /3 in Table|l|are greater than the value (3 = 0.12 
obtained in the classical model with actual qubit poten- 
tial (Sec. IIV B|) . but less than the value j3 — 0.5 in the 
oscillator model [Eq. (|44|l ] . This means that dependence 
on Ti in Eq. (|60l) in the quantum model is in-between 
those found in the oscillator and actual-potential classi- 
cal approaches. Figure 1181 also shows a comparison be- 
tween the results of classical and quantum approaches. 
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FIG. 18: Solid lines: log-log contour plots for the values of the 
error (switching) probability Ps — 0.01, 0.1, and 0.3 on the 
plane of relaxation time Ti (in ns) and coupling capacitance 
Cx (in fF) in the quantum model for (a) N12 = 5 and (b) 
N12 — 10. The corresponding results for Cx,t{Ti) in the clas- 
sical models are shown by the dashed lines (actual potential 
model) and the dotted lines [oscillator model, Eq. 1441 1. The 
numerical data are represented by the points, connected by 
Unes as guides for the eye. The scale at the right corresponds 
to the operation frequency of the two-qubit imaginary-swap' 
quantum gate. 



As one can see, both the oscillator (dotted line) and 
actual-potential (dashed line) classical models give the 
limits for Cx roughly similar to quantum results (within 
a factor of ~ 2), with the exception of the case of the 
oscillator model for iV; = 10 and Ti > 100 ns. The value 
of the coupling capacitance Cx determines the speed of 
two-qubit quantum gates. The right scale of the vertical 
axis in Fig. 1181 converts Cx into the operation frequency 
/i-swap ^ {Cx/C)ujio/2tt of the imaginary-swap£ gate. 

These results can be used for the design of phase-qubit- 
based quantum gates. In particular, they give us the 
maximum allowed coupling capacitance Cx and hence the 
maximum gate operation speed, for a particular tolera- 
ble value Ps of the error due to crosstalk. An important 
result of the quantum treatment is the exponential de- 
pendence of the error probability on Cx and a rather 
slow dependence on Ti. This shows that the measure- 
ment crosstalk is not a big roadblock for the fabrication 
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13 
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0.32 


0.29 


0.28 


0.18 
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TABLE 1: Parameters of Eq. I6OII for the quantum model {Cx 
is in fF and Ti is in ns), limiting the coupling capacitance Cx 
for several values of the error probability Ps and dimensionless 
barrier height Ni2. 



of phase-qubit-based quantum gates with low decoher- 
ence and sufficiently high operation speed. 

The present quantum theory does not take into ac- 
count dissipation in the second qubit. The dissipation 
shortens the effective crosstalk time and thus decreases 
the crosstalk error (the switching probability), similar 
to the classical case discussed in Sec. IIVBI Thus, the 
present results give a lower bound for the maximum al- 
lowed Cx- Taking into account the results of the clas- 
sical model, one may expect ~ 30% larger limit for the 
coupling capacitance (and two-qubit gate frequency) for 
a quantum model with energy dissipation in the second 
qubit. 



VI. CONCLUSIONS 

The main goal of this paper has been to study the 
crosstalk between two capacitively coupled flux-biased 
phase qubits after the measurement pulse. The first 
qubit, which escapes (switches) from the left to the right 
well during the pulse, has been modeled classically. The 
first qubit performs damped oscillations (with energy re- 
laxation time Ti) with increasing in time frequency fd{t)] 
these oscillations perturb the capacitively coupled second 
qubit. The dynamics of the second qubit (which is ini- 
tially in the ground state) has been treated both classi- 
cally and quantum-mechanically. 

In the classical treatment of the second qubit, we have 
compared the previously suggestedS. oscillator model, 
which allows for both analytical and numerical analy- 
sis, with the model based on the exact potential, which 
can be solved only numerically. Both models show a 
sharp resonant excitation of the second qubit. Though 
there is a certain similarity between the two models, they 
significantly differ both quantitatively and qualitatively. 
In contrast to the oscillator model, the exact-potential 
model shows nonlinear and irregular dynamics. The sec- 
ond qubit remains in the left well when the coupling ca- 
pacitance is sufficiently low, Cx < Cx.t, but may escape 
(though not certainly because of complicated dynamics) 
if Cx > Cx,T- We have obtained numerically the de- 
pendence Cx,t{Ti) both in absence and presence of dis- 
sipation in the second qubit for experimentally relevant 
values of the barrier height. 

For the quantum treatment we have developed an effi- 
cient numerical scheme, which uses a subset of eigenstates 
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of the unperturbed Hamiltonian. In this case, similarly 
to the classical case, a fast excitation of the second qubit 
occurs at a moment when the driving frequency fdit) 
is somewhat below the transition frequency between the 
ground and first excited states. However, in contrast to 
the classical case, the switching can now happen even 
when the qubit mean energy is significantly lower than 
the barrier height, either due to tunneling or due to ex- 
citation above the barrier. 

The results for the switching (error) probability Pg 
have been presented as contour plots on the plane of cou- 
pling capacitance Cx and relaxation time Ti (Fig. I18|l . 
Such plots may be important for the design of quantum 
gates based on phase qubits. Comparison of the results 
obtained in the quantum and classical models shows that 
the classical models can be used for a crude estimate 
of the crosstalk error; however, the difference becomes 
significant for Ti ^ 100 ns. In the quantum approach 
the dissipation in the second qubit have been neglected. 



However, by analogy with the classical case, one can ex- 
pect that the account of dissipation will not change the 
results significantly, though it will somewhat increase the 
upper bound for the coupling capacitance, above which 
the crosstalk error becomes intolerable. The model an- 
alyzed in this paper assumes fixed coupling between the 
qubits. Future implementations of experimental schemes 
with adjustable coupling (which can be practically zeroed 
at the time of measurement) will significantly suppress 
the crosstalk error mechanism and correspondingly allow 
for significant increase in the operation frequency of the 
two-qubit quantum gates. 

The present study can also be of relevance for problems 
in other fields (e.g., laser chemistry), which consider ex- 
citation or escape from a potential well by an oscillating 
driving force with changing in time parameters. 

The work was supported by NSA and DTO under ARO 
grant W911NF-04- 1-0204. 
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